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We discuss the relation between bulk topological invariants and the spectrum of surface states 
in three dimensional non-interacting topological insulators. By studying particular models, and 
considering general boundary conditions for the electron wavefunction on the crystal surface, we 
demonstrate that using experimental techniques that probe surface states, only strong topological 
and trivial insulating phases can be distinguished; the latter state being equivalent to a weak 
topological insulator. In a strong topological insulator, only the parity of the number of surface 
states, but not the number itself, is robust against time-reversal invariant boundary perturbations. 
Our results suggest a Z2 definition of the bulk-boundary correspondence, compatible with the Z2 
classification of topological insulators. 
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I. INTRODUCTION 



The defining characteristic of an electric insulator is 
the existence of an energy gap to charge excitations. De- 
pending on the physical origin of that gap insulating ma- 
terials are broadly divided into two classes: Mott insu- 
lators with the gap having an origin in the electron in- 
teractions, and band insulators, where the gap originates 
essentially from the single particle energy terms, with 
many-body effects simply renormalizing the bare band 
parameters. Because of their single-particle nature band 
insulators are often thought to be the simplest systems, 
whose electronic properties are adequately described by 
the usual quantum theory of solids 1 .. Recently, how- 
ever, a classification of these materials has emerged^, 
based on topological invariants 4 , which characterize their 
band structure. In particular, it was shown^ that spin- 
orbit (SO) interaction and time-reversal symmetry can 
stabilize topologically non-trivial electronic states in cer- 
tain systems, which were thus termed topological insula- 
tors (TI). An experimental signature of these phases is 
the presence of chiral metallic surface states, which are 
claimed to be robust against time-reversal invariant lo- 
cal perturbations and weak disorder. Surface states were 
indeed observed in bismuth antimony alloys^, as well as 
in the Bi 2 X 3 and Sb 2 X 3 , with X = (Se, Te), families^. 

The bulk-boundary correspondence is a physical con- 
cept which relates topological properties of the bulk with 
the number of gapless edge modes. This relation has gen- 
erated a lot of discussion in the context of Quantum Hall 
Effect 9 '^ 1 , grapheneii, and now Tfe£ ' 12 i 13 . In particular, 
in Ref. [B| it was argued that a two-dimensional TI is char- 
acterized by an odd number of metallic edge states. Qual- 
itatively, due to time-reversal symmetry, surface states 
in STI do not back-scatter and thus, can not be easily 
localized, becoming robust against local time-reversal in- 
variant perturbations. In three dimensions (3D) it is ar- 
gued that one should distinguish^ among strong and weak 
TI (STI and WTI) states, where the Fermi arc encloses 
an odd or even number of surface Dirac (i.e. Kramers- 
degenerate) points, respectively. 



So far the bulk-boundary correspondence remains a 
quite loosely formulated conjecture. It remains unclear 
how universally valid is this conjecture in real many- 
body insulating systems, or how robust is this correspon- 
dence with respect to variations in the surface proper- 
ties of those systems. In any bulk-boundary correspon- 
dence specific properties of the surface should be of rele- 
vance. Indeed, the electronic edge spectrum can be very 
sensitive 14 to the specific form of the effective surface 
Hamiltonian. For example, in the context of graphene the 
connection between the valley-specific Hall conductivity 
and the number of gapless edge modes has been found to 
be dependent upon the boundary conditions (BC) 1 ^. 

In the present paper we examine the bulk-boundary 
correspondence in 3D band insulators with strong SO 
interaction. We present a technique to analyze their sur- 
face spectra based on the ideas of Ref. Qjl developed in 
the context of semiconductor nanostructures, and study 
the effect of variations in BC on the spectrum of surface 
states. In particular, we use algebraic properties of the 
model Hamiltonian to construct BC according to symme- 
tries of the problem, such as time-reversal invariance. We 
also discuss the role of symmetry-breaking surface pertur- 
bations. Our work leads to a two-fold result. On the one 
hand, we show that even in a clean system with no surface 
reconstruction, from the standpoint of the surface spec- 
tra, there exists no physical distinction between WTI and 
topologically trivial insulating phases. In particular, we 
present examples of trivial insulators, which would ap- 
pear as WTIs, since they also possess an even number of 
surface Dirac cones. On the other hand, the STI state is 
robust against time-reversal invariant boundary pertur- 
bations, in the sense that the Fermi arc always encloses 
an odd number of surface Dirac points with an odd num- 
ber of edge states crossing the Fermi level along a path 
between two time-reversal invariant momenta (TRIMs) 
in the surface Brillouin zone. The number of crossings 
depends on the particular choice of BC. This indicates 
that edge states in a STI are not robust against arbitrary 
time-reversal invariant surface perturbations. However, 
the parity of their number is protected. 

These observations provide a Z 2 formulation of the 
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bulk-boundary correspondence in 3D TI: 

vq = N s mod 2, 

where i/q is the strong Z2 bulk topological invariant^, N s 
is the number of Kramers-degenerate points in the sur- 
face Brillouin zone enclosed by the Fermi arc, and we 
assume that there are no time-reversal breaking pertur- 
bations at the insulator's surface. This relation shows 
that it is only appropriate to differentiate between STI 
and trivial insulating phases, characterized by vq — 1 and 
i>a = 0, respectively. From an experimental point of view 
(by looking at the spectrum of surface states) , a WTI is 
equivalent to a trivial insulator as will be shown. 

The paper is organized as follows. In the next sec- 
tion we consider a simple example of a trivial insulator 
and show that under certain conditions it may possess an 
even number of metallic edge states, thus appearing as 
a WTI. In Sec. Mil we discuss the robustness of surface 
states in a model of a STI, and demonstrate the protec- 
tion of parity of their number, but not the number or 
the states themselves. Moreover, we further discuss the 
indistinguishability between a WTI and a trivial insula- 
tor from the standpoint of the number of edge (surface) 
modes. In Sec. IIV1 we study the effect of time-reversal 
breaking surface perturbations on the edge spectrum of 
a STI. Our conclusions are summarized in Sec. IV] 



II. DIRAC-LIKE BAND MODEL 

The SO interaction plays an important role in es- 
tablishing a band structure characterized by non-trivial 
topological invariants. Physically, this means that proper 
interband matrix elements of the SO term in the Hamil- 
tonian must be comparable with the bulk band gap, so 
that the bands become essentially non-parabolic. As a 
ubiquitous consequence of the band mixing, Dirac-like 
bulk spectra of electrons and holes are formed^. In this 
section we consider a simple model of a topologically triv- 
ial insulator, which neverthless would look like a WTI in 
the sense that the Fermi arc encloses an even number of 
surface Dirac points. We also show that the existence of 
edge states is very sensitive to the BC, imposed on the 
single-particle wavefunction at the crystal surface. 



A. General formalism 

The simplest lattice model, which contains SO interac- 
tion naturally, is just the tight-binding form of the Dirac 
Hamiltonian: 



where from now on subindex i = Xi — [x^, j/,, £j) enumer- 
ates sites of a simple cubic lattice; i + // = //j + 1 with 
|U = x, y, z; = (cicu, Ci va ) destroys an electron with 




FIG. 1. Geometry of the problem studied in this paper. The 
crystal is bounded by the plane z% — and occupies the half- 
space Zi > 0. The unit vector n denotes an outer normal to 
the surface. 



spin a in the conduction (c) or valence (v) band; A and e 
are the SO coupling constant and band gap, respectively; 
and 91 is a U(l) gauge field. The Dirac matrices /3 and 
a M have the well-known properties^: 

{< a v } = 2V,; {of, (3} = 0; [a", a"] = 2ie^ K Z K ; 

with S M = (1 ® <7^) being the 4-spin operator, er^ - the 
Pauli matrix, and - the Levi-Civita symbol. From 
Eq. (TTJ it is easy to obtain an expression for the charge 
current, as a variation of the Hamiltonian with respect 
to the gauge field: 



J? 



5H T 



1 



\Ert a/Pi&j 



(2) 



In the rest of this section, we assume that 9f? = 0, 

The Hamiltonian of Eq. (TTJ commutes with the time- 
reversal (T) and space inversion (V) operators. Also, 
the bulk band structure is invariant under interchange of 
the conduction and valence bands, which means that Hp 
has the charge- conjugation (C) symmetry. These three 
operations are defined by their action on a single-particle 
orbital ip(t, xi) = 4>i(t)^: 

Tip(t, Xi ) =a x a z ip*{-t, Xi ); V^{t, an) = pip(t, -se<); 
Ci>(t,Xi) = pa y ^*{t,Xi). 

Motion in an infinite crystal also preserves the tensor spin 
operato r 16 i 19 : 



r, = ^^(^ 1+fi -h,. 



(3) 



It follows that T is a polar, time-reversal invariant vector. 

In order to study the spectrum of surface states, we 
now consider a half-infinite crystal, bounded by the plane 
Zi = and occupying the region Zi > (see Fig. [T]). On 
the surface the single-particle wavefunction tpf will satisfy 
a linear relation of the form: 



n^^+fj, + Sip is = 0, 



(4) 
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where n is the outer normal to the surface and i s denotes 
lattice sites on the surface. The general structure of the 
4x4 boundary matrix S is fixed by requirements^^ 
that the current (Jo ■ n ) through the surface vanishes, 
and the fundamental symmetries T and V are preserved: 



S = £o/3 + ifi/3a • n. 



(5) 



In this expression the parameters £o.i are free. At 
the phenomenological level, they encode various surface 
properties. For example, the first term in Eq. ([5]) changes 
sign under charge conjugation C, which physically means 
that at the boundary there is a mixing (whose amount is 
controlled by £o) of the bulk Bloch bands. Similarly, the 
second term in S describes the intensity of SO interac- 
tion at the surface, which is caused by rapid changes in 
the crystal field. Thus, after scattering from the surface 
an electron acquires an extra phase, due to spin rotation. 
Since the localized (Tamm) states at the crystal bound- 
ary are formed as a result of interference of bulk Bloch 
states 2 ^, this term has a profound effect on their stability. 
We also note that the BC ([5|) is translationally invariant 
along the surface. 

In the bounded crystal in Fig. [1] only the normal to 
the surface component of T [Eq. ([3])] is conserved. In 
fc-space (fe is the crystal momentum) it can be written 
as: 

T *= *U*i)**(*-i-)**(fci) 



with k± = (k x ,k y ), and 

t z (fcj_ ) = (3 CE X sin k y — Y, y sin k 
The eigenstates of t z (k±) have the form: 

f T k ± (zi)U r k 



(6) 



where 



— ( 1 



t = ±1, s± — y/sin 2 k x -t- sin 2 ky, and the amplitudes tp 
and x are arbitrary. The corresponding eigenvalues are 
ts± . The z-independent part of the kinetic energy in Eq. 
(HJ) can be expressed in terms of t z (k±): 



A 
2i 



h.c. 



= iA 



fe_L ,Zi 



^.Ba z t z {k^ k±z% . 



This fact has two important consequences: (i) The prob- 
lem becomes effectively one- dimensional and we can work 
with two-component wavefunctions, similar to the case of 
a spherically symmetric fielcU^: 



( PTk ± (z l )U Tk± 

XTk ± (z l )U- Tkd 



Xrk ± (z l ) 
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FIG. 2. Stability regions of edge states, obtained from Eq. 
([7]). Inside shaded areas there exist surface states with at 
least one value of the polarization r, as indicated by the num- 
bers ±1. The boundaries of shaded regions are given by the 
expression £i = — (Arsx/e)£o- 



(ii) We are free to choose a representation for the Dirac 
matrices /3 and a z . In the rest of the section, we use the 
following identification: 







B. Surface spectrum 

In order to obtain the edge (surface) spectrum, we fol- 
low the method of Ref. HtJ First, let us consider the bulk 
band structure: 



£ k , ,k. 



± 



A 2 (i 



sm' 



k z ). 



A surface state will have a complex- valued k z = p + iq 
with p and q chosen in such a way that Ime = 0. Since 
q 7^ 0, one obtains: (i) p = or p = ir, and (ii) p — ±ir/2. 
Only in case (i) the state has an energy inside the band 
gap. Therefore, a general (two-component) wavefunc- 
tion, decaying into the crystal (z% > 0), can now be writ- 
ten as a linear combination: 



P— 0,7T 



e + e 



iA(e ip sh<? - ts±) 



Imposing the BC of Eqs. ((4]) and ([5]) one easily obtains 
the inverse localization length q of the state (cf. Ref. Il4h : 



Qh 2 n _ rg T - 1/2 + V(g T - r/2Y + (1 - gj/4 
bn q T k± — 19 i v ' ) 

1 — r z 
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which defines the surface state dispersion relation via 
e Tk± — ±a/£ 2 + ^ 2 ( s ± — sh 2 q T k ± ). In this expression: 



&r(£o,£i) 



<1- 



- - 1 



« + + ! 



G(-l,l). 



A surface state exists only if g > 0, which is equiv- 
alent to g T < 0. Regions in the plane (£o,£i), where 
this inequality is satisfied, are shown in Fig. [2] For a 
given point inside a shaded area, there is a solution ei- 
ther with one, or both values of r. Moreover, there is an 
unshaded region, where the edge spectrum disappears. 
In general, the surface states are gapped. However, this 
gap closes at some points in Fig. [5] For example, let us 
consider the case £o = 0, £i < and r = —rj/^/1 + rf 
with r\ — e/A > 0. Then, solutions for p = and it disen- 
tagle. Therefore, shq T k ± = r), e p Tk ^ = Xc lp rs± = ±Xts±, 
and the wavcfunction is: 



1 

ie ip 



-v - 



T] 2 ] Zl+2 e 



Since £ T k ± depends on k± only through sj_, for the 
chemical potential inside the band gap, the Fermi arc al- 
ways encloses an even number of Dirac points. Thus, 
from an experimental perspective this model system 
would look like a WTI&. Nevertheless, it is straight- 
forward to check that all four topological invariants 
(1/0(^1^2^3)) vanish (we remind the reader that vq = 
and (1/1^2^3) ^ is the mathematical characterization of 
a WTI). As we shall see in the next section, the above 
result is not specific to the model of Eq. (J]) . 



III. LATTICE DIMMOCK MODEL 

The Dirac model ([1} describes a trivial insulator. How- 
ever, it can be extended to support a STI phase. We 
consider one such modification, proposed in Ref. I2H 



H = H 



D 



t 



+ h.c. 



(8) 
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where t > is the hopping amplitude between nearest 
neighbor sites in a 3D cubic lattice. The model thus 
defined is a lattice analog of the well-known Dimmock 
models, which provides an effective-mass description 
of electronic spectra in lead chalcogenides. The Dirac 
model, studied in the previous section, is recovered in 
the limit t — > 0. However, the point t = is singular and 
the limit has to be taken carefully. 

Depending on the ratio of e/t, the model (JSJ exhibits 
the following phases^: (i) STI if 2t < |e| < Gt, (ii) 
WTI for |e| < It and (iii) the trivial band insulator when 
|e| > 6t. In the STI phase the Fermi energy inside the 
gap must cross an odd number of edge states along a path 
between two time-reversal invariant momenta in the sur- 
face Brillouin zone. Indeed, it was shown a long time ago 



that in the long wavelength approximation under band 
inversion the Dimmock model has exactly one Dirac cone 
at the surface^. This conclusion is in agreement with the 
phase diagram, because in the continuum limit the band 
gap A = e — 6t is negative for e < 6t. In this section 
we study the effect of variations in the BC on the sur- 
face spectrum of the lattice Dimmock model, and on the 
physical meaning of the above phase diagram. 

As before, we start by deriving the BC with given sym- 
metry properties. The probability current 



A 



(9) 



vanishes at the surface if the wavefunction satisfies the 
constraint (U) with a boundary operator of the form: 



S 



2t 
A 



io + £,i/3 + i£o/3a • n. 



(10) 



The discussion presented after Eq. ([5]) regarding the 
physical meaning of individual terms is fully applicable 
here as well. 

It is easy to see that the lattice Dimmock model (j8]) 
has the same symmetries as the Dirac model of the pre- 
vious section. Therefore, the problem of determining the 
surface spectrum, in the geometry of Fig. [T] again be- 
comes one-dimensional. However, now it is convenient to 
choose a different representation of the Dirac matrices: 



(11) 



We also introduce a notation e± = e — 2t (cos k x + cos k y ) 
The bulk band structure has the form: 



e k±ikz =±\J (e ± - It cos k z ) 2 + \ 2 (s 2 L + sin 2 k z ) . 

The requirement Ime — is equivalent to the equation 
Ime 2 = under the condition Ree 2 ^ 0. In the com- 
plex fc z -plane (k z = p + iq) Ime 2 vanishes when: (i) 
p = 0, (ii) p — 7r and (iii) along the line cos p dig = 
2te±/[(2t) 2 — A 2 ] . In general, a surface state wavefunc- 
tion is a complicated function with a fcj_-dependent com- 
plex localization length. However, there exists a param- 
eter range, where the edge states are purely evanescent, 
i.e. p = 0. Indeed, let us assume that 

4i<e<10i; {2tf - -(e - ttf ^ A 2 sC (2t) 2 . (12) 
o 

Then, possible values of q are given by: 
2tej_ 



chgi j2 = 



{2t) 2 - A 2 



e\\ 2 + [e 2 - A 2 (si + 1)] [{2t) 2 - A 2 ] 
{2t) 2 - A 2 
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We note that because of peculiar properties of the band 
structure, in the present model there always exists a par- 
tial evanescent solution with chg2 > 1- However, by 
virtue of Eq. (fT2"j) both chgi,2 ^ 1, with qi being the in- 
verse localization length of a state. Using the form (fTTj) 
of a z and /3, the general localized solution can be written 
as: 



E < 

a=l,2 



2ich q a + Ash q a 
e + Xts± 



Before considering generic BC (fit))) , let us analyze the 
case ^i = 0, when there is no band mixing at the bound- 
ary. One can easily show that there exists a surface state 
with energy e T k ± = Xts±, and momentum q 



chgi j2 



2tej 



X 2 e 2 ± -X 2 [(2t) 2 -X 2 



(2t) 2 - A 2 {2tf - A 2 

with the resulting wavefunction 



-qizi 



»»(fo) 



-Q2 Zi 



(13) 



(14) 



where v fa) = [(2t/X-l)^ +e-^]/[(2t/X-l)^+e-^]. 
This state exists only in a region of fc^-space, defined by 
e± < 2t. At the boundary of this region qi = 0, and the 
state merges into the bulk continuum. Clearly, the above 
condition is a priori false when e > 6i, in agreement with 
the general phase diagram 2 ^. 

When £o,i ^ Owe observe that matrices 1, j3 and i/3at-n 
consitute an orthonormal (though not closed) set {m,}: 



Tr mirrij = A5i 



1. Thus, for any value of £i there 



exists a unitary transformation U(uS) = exp (iujot ■ n/2) 
which diagonalizes S (cf. the discussion of the case £i = 
0). For example, in the representation (ITT1) : 



2t 



rtSfa,t 1 )U = -t +i 



X 



Su 



The problem of determining the edge spectrum of the 
original Hamiltonian H [Eq. (|8[l] with the BC Su is 
mathematically identical to the analysis, which led to 
Eq. (IT3l) . Of course, the transformation U changes H as 
well: H — > Hu ^ H, but leaves its topological structure 
intact. Qualitatively, variations in the boundary param- 
eters £o,i can be seen as unitary transformations of the 
bulk Hamiltonian. These transformations also guarantee 
that the BC is kept diagonal. If initially the bulk band 
structure was characterized by a strong topological index 
Mo 7^ 0, the same will be true for Hu- The shape of the 
edge spectrum does depend on a particular choice of £o,i 
and has to be computed explicitly. However, the above 
argument implies that if at £i = 0, the Fermi arc enclosed 
an odd number (e.g. one) of Kramers-degenerate points, 
then the same will be true for any £i- 

In Fig. |31 we present the dispersion relation of sur- 
face states for several values of the BC parameters in 
the STI regime (e = 5.87t) and for the trivial insulator 




FIG. 3. Surface spectra of the lattice Dimmock model for 
different values of the parameters £o,i and e/t. SO inter- 
action is X/t — VS. The spectra are plotted along the path 
(0, 0)-(k x , k y )-(TT, 7r) in the surface Brillouin zone. The shaded 
regions describe the bulk continuum with the band gap equal 
to |e — 6t|. In panel (a) the plot is actually truncated before 
the M-point is reached, for readability. Solid lines represent 
r = +1, dashed lines r = -1. (a) £ = 1-0, £i = 0.0; (b) 
£o = 0.0, Ci = 0.1; (c) Co = 0.0, £i = 0.045; (d) £ = 0.1, 
£i = 0.3. Panels (a)-(c) correspond to a STI with e/t = 5.87, 
panel (d) describes a trivial insulator with e/t = 6.37. 



(e = 6.37i). From panels (a)-(c) it follows that in the STI 
phase the number of edge states at the Fermi level may 
change, depending on the BC parameters, but it always 
remains odd, i.e. the parity of their number is protected. 
These surface states exist for all values of £o,i- 

Panel (d) corresponds to the case e > 6t when the sys- 
tem is a trivial insulator, which can nevertheless exhibit 
metallic edge states. The surface spectrum is quite sen- 
sitive to the BC parameters and disappears for £q.i ~ 1- 
Similar to the problem studied in the previous section, 
the number of Fermi level crossings is even. This circum- 
stance once again shows the physical indistinguishability 
of a WTI and a trivial insulator. 



IV. EFFECT OF TIME-REVERSAL BREAKING 
SURFACE PERTURBATIONS 



Time-reversal symmetry plays a crucial role in stabiliz- 
ing the STI phase and metallic properties of the surface. 
When this symmetry is broken via some physical mech- 
anism, the Kramers theorem no longer holds and surface 
Dirac fermions are expected to acquire mass^^, i.e. they 



become gapped. In Ref. 1241 such gapped edge states were 
indeed observed in Bi2Se3 when doped with magnetic im- 
purities at the surface. 

Phcnomenologically, we can simulate this effect by 
adding a T-breaking perturbation to the boundary oper- 
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ator S. Let us consider one such term: 



AS = irh^oc ■ n, 



(15) 



where 75 = ia x a y a z is the pseudo-scalar matrix-i^. Be- 
cause of this matrix, the correction (fT5|) also breaks space 
inversion symmetry. Nevertheless, it is easy to see that 
the current © vanishes at the surface after AS has been 
added to the boundary operator (TIT)]) . 

Since the operator T z , introduced in Sec. [TTl is invari- 
ant under time-reversal, AS will mix states with different 
values of r. In the absence of a natural conserved quan- 
tity, suitable for labelling single-electron states, the com- 
plete investigation of the edge spectrum becomes quite 
cumbersome. Still, one can understand the qualitative 
effect of AS" by working in the perturbative regime, i.e. 
when ^7- 1 is small. In order to simplify things even fur- 
ther, we confine our analysis to a particular case, when 
there is no band mixing at the surface, i.e. £1 = in 
(fTT)|) . Then, we can use ideas of the previous section to 
transfer the T-breaking term from the BC to the bulk 
Hamiltonian via a unitary transformation U-j-, and em- 
ploy degenerate perturbation theory to treat the result- 
ing corrections in the Hamiltonian within the subspace 
of edge states (fT4|) . 

The unitary transformation lAj- is effected by 75 [cf. 
Eq. GEM: 



U r 



— r i^75/2. 



To first order in £7-, the boundary operator (|10[) remains 
unchanged, and the only correction to the Hamiltonian 
©is: 



AH=^ 
so 



e$> < t r*,-i£(*Jrtt*. /t + h.c. 



1 ,11 



with T = i/3j5. In the basis (fLT| the matrix elements of 
AH can be written as: 



Therefore, at = (Dirac point) the perturbation (fT5|) 
opens a gap A 7- ~ l£r/ £nk in qualitative agreement with 
the experiments of Ref. |2J. 



V. CONCLUSION 

Since in real materials the electronic structure of the 
surface is largerly unknown, investigation of the surface 
properties, such as the spectrum of localized states, re- 
quires a phenomenological approach similar to the one 
employed in the present work. It is important to re- 
alize that any bulk-boundary correspondence cannot be 



simply based on general topological arguments that dis- 
regard hypotheses about the boundary conditions used. 
The key idea of our method consists of using the algebraic 
structure of the bulk Hamiltonian to classify boundary 
conditions for the Bloch wavefunctions according to fun- 
damental symmetries of the problem. 

When the surface preserves time-reversal symmetry, 
our results suggest that experimentally, i.e. by looking 
at the edge spectrum, one can only discriminate between 
the STI and trivial insulating phases. The claimed ex- 
perimental signature of a WTI - an even number of sur- 
face states crossing the Fermi level along a path between 
two TRIMs 5 - turns out to be misleading. Indeed, by 
adjusting phenomenological parameters in the boundary 
conditions, one can make a trivial insulator exhibit a sur- 
face spectrum similar to that of a WTI, even in a clean 
system, where the periodicity of the surface is preserved. 
Therefore, these two phases cannot be physically distin- 
guished by looking at the edge spectrum, and should be 
theoretically classified as the same state. 

On the contrary, in the STI phase the Fermi level 
crosses an odd number of surface states. This precise 
number depends on a particular choice of the boundary 
condition parameters, e.g. it changes from 1 to 3 in pan- 
els (a) and (c) of Fig. El but its parity is always preserved. 
Hence, the surface spectrum of a STI (and for the sake of 
the argument of any band insulator) is not robust against 
time-reversal invariant boundary perturbations. 

These observations can be naturally summarized in a 
consistent definition of a bulk-boundary correspondence: 



N* mod 2, 



where vq is the strong Z2 bulk topological invariant, in- 
troduced by Kane et aZ.— and N s is the number of surface 
Kramers doublets inside the Fermi arc (or the number of 
edge states crossing the Fermi level along any path be- 
tween two TRIMs in the surface Brillouin zone). Only 
the parity of the number of surface states is protected. 
This definition is compatible with the classification of TIs 
according to quantization of the axion #-term&2£, which 
also describes the orbital magneto-electric coupling^. 

The bulk topological invariants are meaningful only if 
the surface preserves time-reversal symmetry. In Sec. IIVI 
we demonstrated that a T-breaking perturbation opens 
a gap in the edge spectrum at the TRIMs, and effectively 
destroys the STI phase. This conclusion is in agreement 
with recent experiments^ performed in bismuth selcnidc 
with the time-reversal symmetry explicitly broken by sur- 
face magnetic impurities. 

Finally, we would like to emphasize that our demon- 
strations are applicable to uncorrelated, i.e. mean-field- 
like, "TP systems. The effect of correlations beyond 
mean-field is still an open problem. While there are some 
efforts^ to extend the existing classification of TIs to in- 
teracting systems, so far they only amount to mean-field 
arguments. 
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